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Recent observations by the Chandra observatory suggest that some neutron stars may cool rapidly, 
perhaps by the direct URCA process which requires a high proton fraction. The proton fraction 
is determined by the nuclear symmetry energy whose density dependence may be constrained by 
measuring the neutron radius of a heavy nucleus, such as '^"^Pb. Such a measurement is necessary 
for a reliable extrapolation of the proton fraction to the higher densities present in the neutron star. 
A large neutron radius in ^°^Pb implies a stiff symmetry energy that grows rapidly with density, 
thereby favoring a high proton fraction and allowing direct URCA cooling. Predictions for the 
neutron radius in '^"^Pb are correlated to the proton fraction in dense matter by using a variety 
of relativistic effective field-theory models. Models that predict a neutron [Rn) minus proton {Hp) 
root-mean-square radius in ^"^Pb to be R„ — Rp < 0.20 fm have proton fractions too small to allow 
the direct URCA cooling of 1.4 AIq neutron stars. Conversely, if Rn — Rp > 0.25 fm, the direct 
URCA process is allowed (by all models) to cool down a 1.4 Mq neutron star. The Parity Radius 
Experiment at Jefferson Laboratory aims to measure the neutron radius in ^"^Pb accurately and 
model independently via parity-violating electron scattering. Such a measurement would greatly 
enhance our ability to either confirm or dismiss the direct URCA cooling of neutron stars. 

PACS numbers: 26.60.-l-c, 21.10.Gv 



I. INTRODUCTION 

Neutron stars are created with very high tempera- 
tures in supernova explosions. Indeed, neutrinos ob- 
served from SN1987A indicate a neutrinosphere temper- 
ature that could be as high as 5 MeV . Neutrons stars 
then cool, primarily by neutrino emission |^ . In the stan- 
dard scenario, the modified URCA reaction, 

n + n ^ n + p + + De, (1) 

emits neutrinos from the volume of the star. This pro- 
cess, however, is relatively slow as a second nucleon is 
necessary to conserve momentum. 

Recent X-ray observations of the neutron star in 
3C58 [||, Vela Q, and Geminga |^ indicate low surface 
temperatures. Moreover, the low quiescent luminosity 
in the transiently accreting binaries KS 1731-260 ||^ and 
Cen X-4 suggest rapid cooling. As X-ray observatories 
progress and our knowledge of neutron-star atmospheres 
and ages improve, additional "cold" neutron stars may 
be discovered. Such low surface temperatures appear 
to require enhanced cooling from reactions that proceed 
faster than the modified URCA process of Eq. (|l|) . En- 
hanced cooling may occur via the weak decay of addi- 
tional hadrons such as pion or kaon condensates 1^, hy- 
perons H, or quark matter (nOl. Yet perhaps the most 
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conservative enhanced-cooling mechanism is the direct 
URCA process [|ll|, ^ of neutron beta decay followed by 
electron capture: 



n 



-P 



(2a) 
(2b) 



This mechanism is not "exotic" as it only requires pro- 
tons, neutrons, and electrons — constituents known to be 
present in dense matter. However, to conserve momen- 
tum in Eq. (|2^) the sum of the Fermi momenta of the 
protons plus that of the electrons must be greater than 
(or equal to) the neutron Fermi momentum. This re- 
quires a relatively large proton fraction. 

Yakovlev and collaborators jl^ are able to reproduce 
measured neutron-star temperatures using a relativistic 
mean-field equation of state that allows direct URCA for 
neutron stars with masses above 1.358 Mq {Mq = so- 
lar mass). In contrast, Tsuruta and collaborators |lj] 
rely on pion condensation to reproduce the measured 
temperatures. They argue that microscopic calculations 
of neutron-rich matter [ p^ using nonrelativistic nucleon- 
nucleon interactions yield too small a proton fraction for 
the URCA process to operate. Unfortunately, these mi- 
croscopic calculations depend on a poorly known three- 
nucleon force and on relativistic effects that could end up 
increasing the proton fraction at high densities. 

Superconductivity and superfluidity can greatly influ- 



ence neutron-star cooling 
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much lower than the pairing gap, pairing correlations 
suppress exponentially the rate of many cooling reac- 
tions. Yet for temperatures of the order of the pairing 



2 



gap, the thermal breaking and subsequent reformation of 
nucleon "Cooper" pairs promotes an additional neutrino- 
emission mechanism that rapidly cools the star [ p^ . 
However, it has been argued in Ref. jl^ that this mech- 
anism alone is unlikely to explain the low temperature of 
some neutron stars. This is because for a large enough 
neutron-pairing gap, pair breaking would rapidly cool all 
neutron stars at a rate almost independent of the mass of 
the star. This would disagree with observations of some 
warm neutron stars. Tsuruta and collaborators have 
claimed that microscopic calculations with a high pro- 
ton concentration show a small proton pairing gap p9| . 
If so, a direct URCA process (one not controlled by pair- 
ing correlations) will cool a star so quickly that thermal 
radiation would become invisible |^|. However, we cau- 
tion that drawing definitive conclusions from microscopic 
calculations of pairing gaps may be premature, as signif- 
icant uncertainties remain in the interactions, equation 
of state, composition, and phases of high-density matter. 

Although the precise mechanism remains unknown, 
some kind of enhanced cooling appears to be required 
to explain the recent observations of cold neutron stars. 
While the need for exotic phases of matter is appealing, 
more conventional cooling scenarios, such as the direct 
URCA process, can not be dismissed on purely theoret- 
ical grounds. Moreover, neutron-star observations alone 
may not be able to resolve the detailed mechanism of en- 
hanced cooling. Thus, we consider complementary lab- 
oratory experiments that could help us confirm (or pos- 
sibly dismiss) the direct URCA process. This can be 
achieved by constraining the symmetry energy of dense 
matter. The symmetry energy describes how the energy 
of (asymmetric) nuclear matter increases as one departs 
from equal numbers of neutrons and protons. The proton 
fraction Yp = Z/A of nuclear matter in beta-equilibrium 
is sensitive to the symmetry energy A large sym- 

metry energy imposes a stiff penalty on the system for 
upsetting the N = Z balance hereby forcing it to retain 
a large proton fraction. 

Energetic heavy-ion collisions probe the symmetry en- 
ergy at high nuclear densities pQ]. Possible observ- 
ables include the ratio of tt^-Io-tt^ production and the 
neutron-proton differential collective flow. However, 
these reactions may suffer from important uncertainties 
associated with the complex strong interactions of the 
heavy-ion collisions. Thus, we rely on a purely elec- 
troweak reaction that can be unambiguously interpreted. 
Parity violating elastic electron scattering from a heavy 
nucleus is sensitive to the neutron density. This is be- 
cause the weak charge of a neutron is much larger than 
the weak charge of a proton. The Parity Radius Exper- 
iment at the Jefferson Laboratory aims to measure the 
neutron radius in ^'^^Pb to a 1% accuracy (±0.05 fm) 
This measurement can be both accurate and model inde- 
pendent . In Ref. Brown showed that the neutron 
radius of *Pb determines the pressure of neutron-rich 
matter at normal densities which, in turn, is related to 
the density dependence of the symmetry energy [M. In 



an earlier work we showed how the neutron radius of 
208pj-, determines properties of the neutron-star surface, 
such as the transition density from a solid crust to a liquid 
interior [^5| . Furthermore, we argued that by comparing 
the neutron radius of ^''^Pb (a low-density observable) to 
the radius of a neutron star (a high and low-density ob- 
servable) evidence may be provided in support of a phase 
transition in dense matter [p6| . 

In the present work we show how the neutron radius 
of a heavy nucleus (such as ^°*Pb) controls the density 
dependence of the symmetry energy. Unfortunately, the 
density dependence of the symmetry energy {dasym/dp) 
is poorly known. Thus a measurement of the neutron 
radius of ^"^Pb seems vital, as it will constrain the den- 
sity dependence of the symmetry energy at low density. 
This, in turn, will allow a more reliable extrapolation 
of the symmetry energy, and thus a more reliable deter- 
mination of the proton fraction at the higher densities 
required in the study of neutron-star structure. While 
in principle collective modes of nuclei, such as the giant- 
dipole or isovector-monopole resonances, are sensitive to 
dagym/dp, in practice this sensitivity is small. Moreover, 
the parameter sets used in the calculations (see various 
tables) have been adjusted so that well known ground- 
state properties remain fixed while changing the neutron 
radius. This shows that existing ground-state informa- 
tion, such as charge densities or binding energies, do 
not determine the neutron radius uniquely. Thus the 
need for a new measurement — such as the neutron ra- 
dius in ^'^^Pb — that will provide important information 
on dasym/dp. 

The paper has been organized as follows. In Sec. ||, 
relativistic effective-field theories for both dense matter 
and finite nuclei are discussed. A large number of param- 
eter sets are considered so that the density-dependence 
of the symmetry energy may be changed while reproduc- 
ing existing ground-state data. In Sec. Ill, results for 



the equilibrium proton fraction as a function of baryon 
density are presented using interactions that predict dif- 
ferent neutron radii in ^"^^Pb. Our summary and con- 
clusions are offered in Sec. In particular, we con- 
clude that for models with a large neutron skin in ^"^Pb 
{Rn — Rp ^ 0.25 fm) the symmetry energy rises rapidly 
with density and the direct URCA cooling of a 1.4 Mq 
neutron star is likely. Conversely, if the neutron radius 
is small {Rn—Rp ^ 0.20 fm) it is unlikely that the direct 
URCA process occurs. In this case, the enhanced cooling 
of neutron stars may indeed require the presence of exotic 
states of matter, such as meson condensates, hyperonic, 
and/or quark matter. 



II. FORMALISM 

Our starting point will be the relativistic effective-field 
theory of Ref. |27| supplemented with additional cou- 
plings between the isoscalar and the isovector mesons. 
This allows us to correlate nuclear observables sensitive 
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to the density dependence of the symmetry energy, such 
as the neutron radius of ^''^Pb, with neutron-star prop- 
erties, such as the threshold mass for URCA coohng. 
As the density dependence of the symmetry energy is 
poorly known, uncertainties in these correlations will be 
explored by considering a wide range of model param- 
eters. The interacting Lagrangian density is thus given 
by 
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(3) 



The model contains an isodoublet nucleon field (V') in- 
teracting via the exchange of two isoscalar mesons, the 
scalar sigma (0) and the vector omega (F^), one isovector 
meson, the rho (b''), and the photon (A'^). In addition 
to meson-nucleon interactions the Lagrangian density in- 
cludes scalar and vector self- interactions. The scalar- 
meson self-interactions (k and A) soften the equation of 
state (EOS) of symmetric nuclear matter at (and near) 
saturation density while the w-meson self-interactions 
(C) soften the high-density component of the EOS. Fi- 
nally, the nonlinear couplings (Ag and Av) are included 
to modify the density dependence of the symmetry en- 
ergy !§. 

The energy of neutron-rich matter may be written in 
terms of the energy of symmetric nuclear matter (pp — Pn) 
and the symmetry energy asym(p)- That is, 

^) = f (P, < = 0) + <'a,y„,(p) + 0{i^) , (4) 
where the neutron excess has been defined as 



t 



Pv 



Pn + Pp 

Here pn is the neutron and pp the proton density, and 



(5) 



P = Pp+ Pn 
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(6) 



The symmetry energy describes how the energy of the 
system increases as one moves away from Pp = Pn- It is 
discussed in Ref. where it is shown that it is given 
by 



asym(/0) 
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where fcp is the Fermi momentum, Ep — \/k'p + M*^, and 
M* =M—gs4>o is the effective nucleon mass. Further, the 
effective rho-meson mass has been defined as follows: 



*2 



ml + 2gl(k,gl<jyl+KglvS) . (8) 



The symmetry energy is given as a sum of two contribu- 
tions. The first term in Eq. (|^) represents the increase 



in the kinetic energy of the system due to the displace- 
ment of the Fermi levels of the two species (neutrons 
and protons). This contribution has been fixed by the 
properties of symmetric nuclear matter as it only de- 
pends on the nucleon effective mass M* . By itself, it 
leads to an unrealistically low value for the symmetry 
energy; for example, at saturation density this contri- 
bution yields ~ 15 MeV, rather than the most realistic 
value of ^ 37 MeV. The second contribution is due to 
the coupling of the rho meson to an isovector-vector cur- 
rent that no longer vanishes in the N ^ Z system. It 
is by adjusting the strength of the NNp coupling con- 
stant that one can now fit the empirical value of the 
symmetry energy at saturation density. However, the 
symmetry energy at saturation is not well constrained 
experimentally. Yet an average of the symmetry energy 
at saturation density and the surface symmetry energy 
is constrained by the binding energy of nuclei. Thus, 
the following prescription is adopted: the value of the 
NNp coupling constant is adjusted so that all parame- 
ter sets have a symmetry energy of asym = 25.67 MeV at 
fci. = 1.15 fm^i(p = 0.10 fm-3) That is. 
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2 (A,g202 ^ AvS^Fo') Afl: 



(9) 
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where Aosym = {asym — kp/(^Ep). This prescription in- 
sures accurate binding energies for heavy nuclei, such as 
^"^Pb. Following this prescription the symmetry energy 
at saturation density is predicted (for Ag = Av = 0) to 
be 37.3, 36.6, and 36.3 MeV for the three families of pa- 
rameter sets considered in this work, namely NL3 [ p8[ , 
S271 ||2^, and Z271[^, respectively (see various tables). 
Moreover, all these parameter sets reproduce the follow- 
ing properties of symmetric nuclear matter: i) nuclear 
saturation at a Fermi momentum of kp = 1.30 fm~^ with 
ii) a, binding energy per nucleon of 16.24 MeV, and in) a 
nuclear incompressibility oi K = 271 MeV. These values 
follow from the successful parametrization of Ref. 
and, thus, have been adopted for the other sets (S271 
and Z271) as well. Yet the various parameter sets differ 
in their values for the effective nucleon mass at satura- 
tion density, ii) the value of the w-meson quartic coupling 
(C 7^ for set Z271 but vanishes for the NL3 and S271 
sets) and Hi) the nonlinear couplings Ag and Av (see var- 
ious tables). Note that changing Ag or Av modifies the 
density dependence of the symmetry energy through a 
change in the effective rho-meson mass [see Eq. (^)]. In 
general, increasing either Ag or Av causes the symmetry 
energy to soften, that is, to grow slower with increasing 
density. This, in turn, allows for a larger neutron-proton 
mismatch or equivalently, for a lower equilibrium proton 
fraction at high density. The neutron radius of ^'^^Ph 
also depends on the density dependence of the symme- 
try energy. A stiff density dependence {i.e., pressure) 
for neutron matter pushes neutrons out against surface 
tension, leading to a large neutron radius. The pressure 
of neutron-rich matter depends on the derivative of the 
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TABLE I: Model parameters used in the calculations. The 
parameter k and the scalar mass rris are given in MeV. The 
nucleon, rho, and omega masses are kept fixed at M = 939, 
nil, — 763, and = 783 MeV, respectively — except in the 
case of the NL3 model where it is fixed at rriuj = 782.5 MeV. 



Model 




a! 






A 


c 


NL3 


508.194 


104.3871 


165.5854 


3.8599 


-0.01591 


0.00 


S271 


505.000 


81.1071 


116.7655 


6.6834 


-0.01580 


0.00 


Z271 


465.000 


49.4401 


70.6689 


6.1696 


+0.15634 


0.06 



TABLE IL Results for the NL3 parameter set with As = 0. 
The neutron skin {R„ — Rp) of '^"^Pb is given along with 
the threshold density for the direct URCA process Purca, 
the corresponding proton fraction i^uRCA, and the minimum 
mass neutron star where the direct URCA process is allowed 
A^uRCA- Note that Rn — Rp is given in fm, Purca in fm~^, 
and MuRCA in solar masses. 



Av 


9l 


Rn—Rp 


Purca 


VpURCA 


A/uRCA 


0.0000 


79.6 


0.280 


0.205 


0.130 


0.838 


0.0050 


84.9 


0.266 


0.233 


0.131 


0.944 


0.0100 


90.9 


0.251 


0.271 


0.132 


1.224 


0.0125 


94.2 


0.244 


0.293 


0.133 


1.435 


0.0150 


97.9 


0.237 


0.319 


0.134 


1.671 


0.0200 


106.0 


0.223 


0.376 


0.135 


2.123 


0.0250 


115.6 


0.209 


0.442 


0.136 


2.449 



energy of symmetric matter with respect to the density 
{dE{p,t = 0)/dp) and on the derivative of the symme- 
try energy (dasym/dp). While the former is well known, 
at least in the vicinity of the saturation density, the lat- 
ter is not. Hence, by changing the values of Ag or Av 
one can adjust the density dependence of the symmetry 
energy dagym/dp, while keeping a variety of well-known 
ground-state properties (such as the proton radius and 
the binding energy of ^'^^Pb) unchanged. Note that pa- 
rameter sets with a large "pressure" , dasym/dp, yield a 
large neutron radius in ^°*Pb. 



III. RESULTS 

In this section results are presented for various observ- 
ables that have been computed using an equation of state 
for matter composed of neutrons, protons, electrons, and 
muons in beta equilibrium: 



n 



p + e H- i^e , 

^ p,^ + l^e + V^i 



(10a) 
(10b) 



These reactions demand that the chemical potential (/i) 
of the various constituents be related by the following 
equation: 



TABLE IIL Results for the S271 parameter set with As = 0. 
The neutron skin (_R„ — Rp) of ^"^Pb is given along with 
the threshold density for the direct URCA process Purca, 
the corresponding proton fraction I^urca, and the minimum 
mass neutron star where the direct URCA process is allowed 
AfuRCA- Note that Rn—Rp is given in fm, Purca in fm~^, 
and MuRCA in solar masses. 



Av 


9l 


Rn — Rp 


Purca 


Yp URCA 


AfuRCA 


0.000 


85.4357 


0.254 


0.224 


0.130 


0.830 


0.005 


88.3668 


0.246 


0.252 


0.132 


0.894 


0.010 


91.5061 


0.238 


0.296 


0.133 


1.059 


0.015 


94.8767 


0.230 


0.374 


0.135 


1.429 


0.020 


98.5051 


0.221 


0.501 


0.137 


1.938 


0.025 


102.4221 


0.214 


0.663 


0.139 


2.248 


0.030 


106.6635 


0.205 


0.843 


0.140 


2.343 



TABLE IV; Resufis for the Z271 parameter set with As = 0. 
The neutron skin [Rn — Rp) of '^"^Pb is given along with 
the threshold density for the direct URCA process Purca, 
the corresponding proton fraction I^urca, and the minimum 
mass neutron star where the direct URCA process is allowed 
AfuRCA- Note that Rn—Rp is given in fm, Purca in fm~^, 
and A/uRCA in solar masses. 



Av 


9l 


Rn — Rp 


Purca 


Yp URCA 


AfuRCA 


0.000 


90.2110 


0.241 


0.242 


0.131 


0.816 


0.010 


92.5415 


0.235 


0.274 


0.132 


0.862 


0.020 


94.9956 


0.228 


0.332 


0.134 


0.971 


0.025 


96.2721 


0.225 


0.386 


0.135 


1.079 


0.030 


97.5834 


0.222 


0.500 


0.137 


1.270 


0.035 


98.9310 


0.219 


0.747 


0.139 


1.498 


0.040 


100.3162 


0.215 


1.028 


0.141 


1.583 



where pm fJ-p^ Me, and p^ represent the chemical poten- 
tials of neutrons, protons, electrons, and muons, re spec - 
tively. Neglecting the rest mass of the electron, Eq. (10b) 
is equivalent (for fcp >mp) to the following equation ex- 
pressed in terms of the Fermi momenta of the electron 
and the muon: 



(12) 



Finally, charge neutrality imposes the following con- 
straint on the system: pp = pe + Pfj., or equivalently. 



(13) 



(11) 



For a given proton Fermi momentum, the correspond- 
ing Fermi momenta for the electrons and the muons are 
readily obtained by solving Eqs. ( p^ ) and (13). With 
these in hand, Eq. ( pj]) determines the equilibrium neu- 
tron {Yn = N/A) and proton {Yp = Z/A) fractions in the 
system. 

In Fig. |l| the proton fraction Yp for matter in beta 
equilibrium is shown as a function of the baryon den- 
sity for all models discussed in the text (see Table |). 
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FIG. 1: Proton fraction Yp versus baryon density p for neutron-rich matter in beta equilibrium for all parameter sets discussed 
in the text. For every given set, the curves represent different values of the isoscalar-isovector coupling (Av or As) that yield 
the values of 7?„ 



indicated in the inset. See also the various tables in the text. 



TABLE V: Results for the Z271 parameter set with Av = 0. 
The neutron skin (_R„ — _Rp) of ^"^Pb is given along with 
the threshold density for the direct URCA process Purca, 
the corresponding proton fraction VpURCA, and the minimum 
mass neutron star where the direct URCA process is allowed 
A^uRCA- Note that Rn — Rp is given in fm, Purca in fm~^, 
and MuRCA in solar masses. 



As 


9l 


Rn—Rp 


Purca 


Yp URCA 


A^URCA 


0.000 


90.2110 


0.241 


0.242 


0.131 


0.816 


0.010 


96.3974 


0.229 


0.287 


0.133 


0.901 


0.020 


103.4949 


0.216 


0.366 


0.135 


1.078 


0.030 


111.7205 


0.204 


0.488 


0.137 


1.300 


0.040 


121.3666 


0.191 


0.636 


0.139 


1.467 


0.050 


132.8358 


0.178 


0.789 


0.140 


1.560 


0.060 


146.6988 


0.164 


0.936 


0.141 


1.605 



The various curves displayed in each panel are for val- 
ues of As or Av that give the indicated values for the 
neutron skin of ^°*Pb. Note that the neutron skin of a 
nucleus is defined as the difference between the neutron 



(Rn) and the proton (Rp) root-mean-square radii. All of 
the curves yield the same proton fraction at low density 
(fcp = 1.15 fm~^ or p « 0.1 fm~'^) because the symme- 
try energy has been adjusted to agym = 25.67 MeV in 
order to reproduce the binding energy ^°*Pb. The pro- 
ton fraction increases more rapidly with density for those 
parameter sets that yield larger neutron radii in ^ospj-, 
(namely, those with a stiffer symmetry energy). Thus, 
the neutron radius of ^"^Pb constrains the slope dYp/dp 
at normal densities. This enables one to make a more 
reliable extrapolation of Yp to the higher densities where 
it displays some model dependence. 

The direct URCA process is viable only when the pro- 
ton fraction is large enough to conserve momentum in 
the neutron beta decay reaction of Eq. (|2a|). Hence, for 
this reaction to proceed, the Fermi momenta of neutrons, 
protons, and electrons must satisfy the following relation: 

k^<kP+k^. (14) 

The URCA threshold density Purca is defined as the den- 
sity at which the equality {kp = kp+kp) is satisfied. Note 
that in the simplified case of matter without muons, that 
is, fcp < and kp = k^, the proton fraction at the onset 
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of the direct URCA process is = 1/9 ~ 0.111. In the 
opposite hmit, fcp^m^, the threshold proton fraction is 
l^'^ 0.148. Thus, the threshold proton fraction must be 
contained within these two values for all baryon densities 
(see Tables In Fig. | the URCA threshold density 

is displayed as a function of the neutron skin i?„ — Rp of 
208pb. There is a clear tendency for Purca to decrease 
with increasing neutron skin. Recall that a large neutron 
skin implies a stiff symmetry energy and a large proton 
fraction. Thus the onset for the direct URCA process 
for models with large neutron skins occurs at low baryon 
densities. Indeed, parameter sets with neutron skins of 
Rn — Rp ^ 0.24 fm yield a relative low URCA density 
of PuRCA ^ 0.30 fm . This density is only a factor of 
two larger than normal nuclear matter saturation density 
(po ~ 0.15 fm~^). In contrast, if Rn — Rp < 0.21 fm, the 
onset for the URCA process is above 3po- 
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in models with large neutron skins occurs at low baryon 
densities, thus lower "URCA masses". If the neutron 
skin of 208p]3 jg ^ggg tj^ai^ about Rn~Rp < 0.20 fm, then 
all parameter sets considered in this work predict that 
a neutron star of 1.4 Mq will not undergo URCA cool- 
ing. Conversely, if Rn—Rp > 0.25 fm, then all parameter 
sets allow URCA cooling for 1.4 Mq neutron star. Note 
that all well measured neutron stars have masses near 
1.4 Mq The threshold neutron star mass for the di- 
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FIG. 3: Threshold neutron-star mass for the direct URCA 
process versus the predicted neutron skin (Rn-Rp) of ^"^Pb. 
Parameter sets NL3, S271, and Z271v use a nonzero value for 
Av while Z271s uses a nonzero As. 



FIG. 2: Threshold density for the direct URCA process versus 
the predicted neutron skin {R„-Rp) of ^"^Pb. Parameter sets 
NL3, S271, and Z271v use a nonzero value for Av while Z271s 
uses a nonzero Ag. 

The structure of spherical neutron stars in hydrostatic 
equilibrium is solely determined by the equation of state 
of neutron-rich matter in beta equilibrium. Having spec- 
ified the equation of state, we determine the mass of neu- 
tron stars that may cool via the direct URCA process by 
integrating the Tolman-Oppenheimer-Volkoff equations. 
Our treatment of the low-density crust, where the mat- 
ter in the star is nonuniform, will be discussed in a later 
work This region, however, has almost no effect 

on neutron-star masses. As mentioned earlier, we con- 
sider matter composed of neutrons, protons, muons, and 
(massless) electrons. In Fig. ^ we display, as a function 
of the neutron skin in ^°*Pb, the mass of a neutron star 
whose central density equals the URCA density (purca) 
of Fig. ^. Neutron stars with larger masses, and thus 
higher central densities, will cool by the direct URCA 
process; those with lower masses will not. There is an 
obvious trend for this threshold mass to decrease with in- 
creasing Rn—Rp- Recall that the onset for URCA cooling 



rect URCA process, Murca, depends on both the "criti- 
cal" URCA density Purca (of Fig. ||) and on the equation 
of state at low and high densities. Yet the neutron skin of 
^•^^Pb constrains only the low density EOS j23|, thereby 
generating the model dependence displayed by Fig. ^. In 
particular, the radius of a neutron star, although corre- 
lated to Rn — Rp, is not determined uniquely by it [ p6[ . 
This suggests that models with a stiff EOS at high den- 
sity, such as the NL3 parameter set, will yield neutron 
stars with relatively large radii and low central densities. 
This implies, for example, that for a fixed central density 
of PuRCA = 0.3 fm~'^, the NL3 set (with large radii) gen- 
erates an URCA mass of Murca — 1.4 Mq; in contrast, 
the softer Z271s set (with small radii) yields an URCA 
mass of only Murca — 1 Mq. Equivalently, to make a 
Murca = 1.4 Mq neutron star the NL3 set requires an in- 
terior density of Purca = 0.3 fm~"^, while a central density 
almost twice as large is needed for the Z271s set to gener- 
ate the same mass neutron star. These facts suggest that 
an accurate measurement of neutron-star radii may re- 
duce the model-dependence (i.e., the spread) observed in 
Fig. I Yet, even without further constraints the spread 
is relatively small and a measurement of Rn—Rp in ^"^Pb 
may still prove decisive. 
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IV. SUMMARY AND CONCLUSIONS 

Recent X-ray observations suggest that some neutron 
stars cool quickly. This enhanced cooling could arise from 
the direct URCA process — that requires a high proton 
fraction — or from the beta decay of additional hadrons in 
dense matter, such as pious, kaons, hyperons, or quarks. 
Yet, it seems unlikely that the X-ray observations alone 
will determine the origin of the enhanced cooling. 

In this work we propose to use a laboratory experiment 
to constrain the direct URCA process in neutron stars. 
The Parity Radius Experiment at the Jefferson Labora- 
tory [ pT| , P4 aims to measure the neutron radius of ^"^Pb 
accurately and model independently via parity-violating 
electron scattering. For the direct URCA process to be 
realized, the equilibrium proton fraction in the star must 
be large. The equilibrium proton fraction is determined 
by the symmetry energy, whose density dependence can 
be strongly constrained through a measurement of the 
neutron radius in ^"^Pb. Such a measurement could pro- 
vide a reliable extrapolation of the proton fraction to 
higher densities. Thus, predictions for the neutron radius 
in ^"^Pb have been correlated to the proton fraction in 
dense neutron rich matter by using a wide range of rela- 
tivistic effective-field theory models. We find that models 
with a neutron skin in ^'^^Pb of Rn—Rp < 0.20 fm gener- 
ate proton fractions that are too small to allow the direct 
URCA process in 1.4 Mq neutron stars. Conversely, if 
Rn — Rp ^ 0.25 fm, then all models predict the URCA 
cooling of 1.4 Mq stars. 

While this paper has focused on rclativistic effective 
field-theory models, we expect our conclusions to be gen- 
eral and applicable to other approaches, both rclativis- 
tic and nonrelativistic. For example, the nonrelativistic 
equation of state of Friedman and Pandharipande 



predicts too small a proton fraction for URCA cooling 
to be possible. Moreover, this equation of state yields a 
neutron skin in ^ospb of only i?„-i?p = 0.16±0.02 fm [||. 
Thus, these results are fully consistent with Fig. ^ that 
predicts no URCA cooling for such a small value of 
R71 Rp • 

The equation of state considered in this work consists 
of matter composed of neutrons, protons, electrons, and 
muons in beta equilibrium; no exotic component was in- 
voked. Further, no explicit proton or neutron pairing was 
considered. Nucleon superfluidity is an accepted phe- 
nomenon in nuclear physics and superfluid gaps are im- 
portant for the cooling of neutron stars ||l^. Thus, the 
study of pairing gaps in relativistic effective-field theo- 
ries is an important area of future work; first steps in 
this direction have been taken in Ref. ||31| . In particu- 
lar, the proton pairing gap in matter with a high proton 
concentration must be computed jl^ . 

In summary, the feasibility of enhanced cooling of neu- 
tron stars via the direct URCA process was studied by 
correlating the proton fraction in dense, neutron-rich 
matter to the neutron skin of ^°^Pb. Thus, a measure- 
ment of the neutron radius in ^"^^Pb may become vital for 
confirming (or dismissing) the direct URCA cooling of 
neutron stars. If direct URCA cooling is ruled out, then 
observations of enhanced cooling may provide strong ev- 
idence in support of exotic states of matter, such as me- 
son condensates and quark matter, at the core of neutron 
stars. 
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